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Abstract 



^ . The transverse broadening of an energetic jet passing through a non-AbeUan plasma is be- 

_j. . heved to be described by the thermal expectation value of a light-cone Wilson loop. In this 

■^ , exploratory study, we measure the light-cone Wilson loop with classical lattice gauge theory 

• ! simulations. We observe, as suggested by previous studies, that there are strong interac- 

f~^ \ tions already at short transverse distances, which may lead to more efficient jet quenching 

^^ \ than in leading-order perturbation theory. We also verify that the asymptotics of the Wilson 

loop do not change qualitatively when crossing the light cone, which supports arguments in 
J ^ the literature that infrared contributions to jet quenching can be studied with dimension- 

'V^ I ally reduced simulations in the space-like domain. Finally we speculate on possibilities for 

.?.' full four-dimensional lattice studies of the same observable, perhaps by employing shifted 

boundary conditions in order to simulate ensembles boosted by an imaginary velocity. 
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1. Introduction 

When an energetic jet traverses a strongly interacting thermal medium, various interactions 
take place and lead to dissipation: the jet loses some of its energy and sharpness. The latter 
phenomenon is referred to as jet broadening, or jet quenching. If its efficiency is measured 
experimentally as a function of the jet's energy (this can be done particularly well if the total 
jet momentum is balanced against that of a hard photon, which does not lose energy to the 
medium [1]), then we may learn something about the properties of the medium itself. The 
current understanding is that in order to explain the jet quenching observed empirically in 
heavy ion collision experiments, interactions have to be much stronger than suggested by 
leading-order perturbation theory (for reviews see, e.g., refs. [2]-[8]). 

On an intuitive level, a highly energetic jet can be thought of as a light-cone Wilson line, 
and the fact that we are probing its fate in the transverse direction leads us to correlate the 
Wilson line with a slightly displaced Hermitean conjugate. Adding lines at both ends leads 
to a light-cone Wilson loop. Arguments have been given to make the correspondence precise 
(see, e.g., refs. [3],[9]-[12]), however it appears difficult to state the form of the error that is 
made in this approximation. In the following we take the light-cone Wilson loop as a starting 
point, without dwelling any further on its relation to physically measurable quantities. 

In a statistical environment (with a temperature T, assumed to be above a few hundred 
MeV), thermal noise leads to decoherence. As a result the light-cone Wilson loop, to be 
denoted by W, "decays" at large Minkowskian times t » h/T} Schematically, assuming an 
appropriate time ordering, we may expect that 

{W{t,r^)) '^^/"^ Z(rx)e-'^(^^)* ~ ^(^^) e-Rey(r,)i g-limy(r,)|i ^ (^^1) 

where r±_ = \t±_\ is the length of a 2-dimensional transverse vector; Rey(r^) is a real phase; 
and {...)rp refers to a thermal expectation value. If the coefficient of the exponential decay is 
represented in Fourier space, 



|ImF(rx)| = / {l-e^^^-^^)C{k^), (1.2) 

then C{k±) is often referred to as the "transverse collision kernel" ([13] and references therein). 
Considering for concreteness a Wilson loop in the fundamental representation, the leading- 
order expression for C{k±) at small transverse momenta reads ([14], eq. (44)) 

Cik,,^iTC,[''yo[?^). (1.3) 



^ Since the concept of a classical limit appears frequently, it is useful to show h explicitly, thereby keeping 
the units of time and energy separate. In contrast we set the speed of light equal to unity as usual. 



where g^ = Anas/h is the strong gauge couphng; Cp = {N^ — l)/(2iVc); and 

is the Debye mass parameter (which has units of inverse distance squared). The question 
we are interested in is how large the corrections to eq. (1.3) can be, particularly within the 
infrared domain k± <^ ttT/H. 

Previous work already exists on infrared corrections to eq. (1.3). In particular, the cor- 
rections of 0{g'^T'^) were computed for k±_ ~ m^ in ref. [13], and non-perturbative effects 
of 0{g^T'^) for k^ ~ Q^T/ir were addressed in ref. [15]. In ref. [13] it was noted that for 
k± ~ TTT-E the perturbative series might be slowly convergent, and therefore in need of an 
all-orders resummation. Conceptually, the aim of the current study is to implement such a 
resummation through numerical simulations of a low-energy description. 

More precisely, our goal is to address eq. (1.1) within the framework of classical lattice 
gauge theory (CLGT). It should be immediately acknowledged that although CLGT does 
represent^ the physics of the system at scales A;j_ ~ g'^T/ir, it actually is not quantitatively 
accurate at the scales k± ~ m^ that are of most interest here. The reason is that it is 
highly sensitive to lattice artifacts in this momentum range [23, 24]. Nevertheless, it still 
contains the correct physics on the qualitative level; indeed CLGT simulations have been 
useful for gaining insight on various phenomena at the Debye scale (see e.g. refs. [25, 26]), 
thereby serving as a stepping stone towards full four-dimensional simulations of the same 
problems (see e.g. refs. [27]-[31]). The great strength of CLGT is that it operates directly in 
Minkowskian space-time, thereby circumventing all issues related to analytic continuation. 
The purpose of the present study is to explore what CLGT can teach us about the light-cone 
Wilson loop in the domain indicated in eq. (1.1).^ 

The plan of this paper is the following. After outlining the general framework (sec. 2), we 
present some analytic expectations in sec. 3, setting the stage for a comparison with numerical 
data. The numerical results are presented in sec. 4, and we conclude in sec. 5. 



2. General framework 

With a view on obtaining a formulation which may eventually be amenable to full four- 
dimensional lattice Monte Carlo simulations, we start by defining a "tilted" Wilson loop 



■^Originally CLGT simulations were employed for addressing the rate of non-perturbative anomalous chi- 
rality violation originating from the scale k±_ ~ g^T /-k, see e.g. refs. [16]-[21] and references therein. They 
have also been used for studying the dynamics of thermal phase transitions, see e.g. ref. [22], as well as many 
non-equilibrium problems in cosmology and heavy ion collision experiments. 

^Previously CLGT simulations have been used as an ingredient in a phenomenological study of jet quenching 
of hard particles [32], but the light-cone Wilson loop was not measured. 
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Figure 1: An illustration of a tilted Wilson loop in Euclidean space-time. The slope in the "parallel" 
direction (r||) is parametrized by a Euclidean velocity v-e, so that the right edge lies at r|| = v^t. The 
transverse extent {r±) can be interpreted as the length of a two-dimensional vector. 



in Euclidean space-time. The Wilson loop is parametrized by a transverse extent, rj_; by 
an imaginary-time variable, r G (0, /3), where /3 = h/T; and by a velocity, v^- At the 
end of the computation both r and v^ will be subjected to a Wick rotation, but for the 
moment they are treated as real variables. The Wilson loop is illustrated in fig. 1. In the 
limit Ve — ?• 0, it goes over into the Wilson loop defined in the context of heavy quarkonium 
physics in ref. [33]. (We note that it may ultimately be more useful to "tilt" the thermal 
ensemble rather than the Wilson loop, which in four dimensions can be achieved through 
shifted boundary conditions [34].) 

More concretely, starting with the continuum formulation and choosing sign conventions 
in which the covariant derivative in the fundamental representation is D^ = d^ + ig^A^, a 
straight Wilson line reads 



W[X2;X^]=VeM-m [ 'dX^A^), (2.1) 

JXi 

where X = (r, x). The foremost tilted line of fig. 1 can be expressed as 

VF[(r,rx+VEr);(0,rx)] = 1 - i^o / dn (Aq + v^ • A)(Ti,r^ + VeTi) + . . . , (2.2) 

Jo 

where Vg = v^e^^. The expectation value of the Wilson loop is defined as 

CE(T,VE,r^) = ^Tr(w[{0,Oy, {T,^rET)■, (r,r^+VEr); (0,r^); (0,0)])^, (2.3) 

where the thermal average {■■■)rp implies periodic boundary conditions for bosonic and an- 
tiperiodic ones for fermionic fields over the Euclidean time direction. In the following, we 
have in mind evaluating the expectation value within pure SU(3) gauge theory, even though 
this restriction can in principle be relaxed. 

Since the physical observable that we are interested in refers to Minkowskian time, an 
analytic continuation needs to be carried out at the end of the computation. Technically, 



we do this by substituting r — )■ it, which for 2-point functions yields the time ordering 
corresponding to a Wightman correlator denoted by C>. (The Wilson loop can always be 
thought of as a 2-point function in time if the tilted lines are gauged to unity; general issues 
related to time ordering have been discussed in refs. [3],[9]-[13].) However, since in the 
following we will simultaneously take the classical limit, time ordering actually plays no role. 
The classical limit can be defined by writing 

/3 = ^ , 9i= g'h, (2.4) 

and subsequently setting ?i — )■ [35] . This limit is non-trivial and results in an interacting non- 
Abelian gauge theory which captures the infrared features of the system's real-time thermal 
dynamics [16, 17]. 

Apart from the continuum formulation, we also consider a lattice formulation of the theory 
in the following. Like in ref. [36], the theory is discretized only in spatial directions, with a 
finite lattice spacing a, whereas the time direction remains continuous. Thereby the four- 
dimensional Euclidean action can formally be expressed as 

^E ^ a^j; rdr{X;Tr[i?,2(X)] + ^X;TY[l-P,,(X)]}, (2.5) 

where gQ denotes the bare gauge coupling and Ei, Pij denote the electric field strength and 
the spatial plaquette, respectively: 

^^(^) ^ ^[^rU,{X)]uj{X) ^ AoiX)-U^{X)Ao{X + ae,)uJ{X) ^ ^^^^^ 

* a^o a ' 

P^J{X) = U^{X)Uj{X + aei)Ul{X + ae,)u]{X). (2.7) 

Here U,i G SU(3) are link matrices, and ^o is a traceless and Hermitean gauge field. The 
action is invariant under the gauge transformation 

Ui{X) -^ G{X)Ui{X)G^\X + aei) , (2.8) 

A^{X) ^ G{X)A^{X)G-\X) + -[drG{X)]G-\X), (2.9) 

% 

with G G SU(3). For perturbative computations we make use of covariant gauges; in contrast, 
on the real-time simulation side it is convenient to make use of the corresponding Hamiltonian 
formulation with a vanishing Minkowskian Aq and a corresponding Gauss law constraint. 



^This formulation is invoked because of its close relation to CLGT; in contrast, the speculations to be 
made about full four-dimensional lattice studies in sec. 5 apply equally well to the standard formulation with 
a symmetric discretization in all directions. 
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Figure 2: The graphs contributing to the tilted Wilson loop of fig. 1 at ©(ffo)- Wiggly lines stand 
for HTL-resummed gluon propagators. 

3. Analytic expectations 

3.1. HTL result in continuum 

Our ultimate goal is to compute the analytic continuation of eq. (2.3) at large Minkowskian 
times, t > TT/{g'^T), and large transverse distances, r± > l/m-E- We start, however, by inspect- 
ing short distances, r^<l/m^. This can be done with perturbation theory, provided that 
we recall that at high temperatures the loop expansion needs to be resummed to all orders 
in order to arrive at a consistent weak-coupling result. We are working at leading non-trivial 
order in this regime, and then the effects of resummation are contained within Hard Thermal 
Loop (HTL) [37, 38] propagators. 

Concretely, we carry out the computation by evaluating the graphs of fig. 2 with the 
Euclidean propagator 



Al{X)Al{Y)) = d^^i e^^(^-^) 



K 



Wl,{K) ^ Wl,{K) ^ iK,K, 



i^2 + n K2 + n 



K^ 



(3.1) 



where K = (A;„,k) and ^ is a gauge parameter. The projectors read 






(K) 



s^Aj [ ^ij 



fxi'i rvj 

1^ 



n^iK) = 5^ 



flU 



K2 






{K) 



(3.2) 



The Euclidean propagators are expressed in a spectral representation, 

1 /■°°dA:oPT(E)(^) 



K2 + H 



t(e) 



vr kn 



LlXiy^ 



(3.3) 



where /C = (feo, k), and subsequently the Matsubara sums are carried out. The explicit forms 
of the self-energies can be found in the literature but are not needed here. In general the 
computation parallels that in ref. [33], except that it is in some sense simpler (as long as 
we stay in continuum) : indeed a non-zero f e "regulates" the contributions of the Matsubara 
zero modes, so that they no longer need to be treated separately from the non-zero ones. 



Some technical details of the computation are presented in appendix A. Here we merely 
note that analytic continuation is carried out as r — t- it, v^ — ?• —iv, and the classical limit is 
taken as discussed around eq. (2.4).^ The definition of a potential reads (cf. eq. (1-1)) 

idtCj^{it, -iv, r_i_) = V{t, v, r^)C^{it, -iv, r^) ; (3.4) 

taking the limit t — )■ oo and setting f — t- 1, we reproduce the result of eq. (1.3): 



Kr^(oo,l,rx) = -.<7^TCpy^^(l-cosk^.r^)^^-^^-^-^) (3.5) 



-I 

2-K 



Inf ^ ^ ) + 7e + KQ{m^ri_] 



(3.6) 



Here Kq is a modified Bessel function. 

Next-to-leading order (NLO) corrections to the integrand of eq. (3.5) have been determined 
in ref. [13]. They are large and increase the magnitude of the imaginary part; their numerical 
contribution to eq. (3.6) is shown in fig. 4 below. 

3.2. HTL result on a spatial lattice 

For a practical measurement, the theory needs to be regularized; within CLGT, this means 
that we consider (a Minkowski-space classical limit of) the theory defined by eq. (2.5). Ex- 
pressing everything in lattice units and taking the limit of eq. (2.4), the results depend on a 
single parameter, which we denote by 

/3g - 4^ . (3.7) 

Initial configurations are generated with the weight exp(— /^c-ffci) nx'^(^(-'^))' where 

3 , 3 






(3.8) 



^(x) denotes the Gauss law constraint; and £i{x.) are suitably normalized canonical momenta 
conjugate to the link matrices Ui{x.). Subsequently the fields are evolved according to classical 
equations of motion (cf. eqs. (4.1), (4.2)), and the observable is measured as illustrated in 
fig. 3. (Further details on CLGT simulations can be found e.g. in refs. [17]-[20], [25, 26]; the 
normalization of the electric field is strongly reference-dependent.) 

Within CLGT, the Debye mass scale of the continuum formulation gets replaced with 
7TT-I — 7- ff^T /a, whereas the coupling constant scale remains put at g^T. In lattice units, this 



^In practice the classical limit amounts to assuming that fco <^ -kT /h; therefore, at leading order it correctly 
represents the physics of the large-time or low-energy limit of the exponential decay. 
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Figure 3: An illustration of a tilted Wilson loop after the discretization of the spatial directions with 
a lattice spacing a, and a Wick rotation of both t and v-^ to Minkowskian space-time. We represent 
the tilted Wilson lines by averaging the smallest possible building blocks over the upper and lower 
paths. Measurements are taken at values t — na/v, with n € N. 

implies that we want to determine the Wilson loop at separations r±/a> \//3g and time 
scales t/a>l3G- The latter of these requirements poses a significant challenge at large /3g, 
and introduces a source of systematic errors with any limited resources. Approaching this 
regime from below, perturbation theory can again be used, but necessitates a HTL-type 
resummation, whose details were worked out in refs. [23, 24]. 

In practice, carrying out perturbative computations even to leading non-trivial order is 
cumbersome, due to the asymmetry in the discretizations of the temporal and spatial direc- 
tions. As an example, the expression obtained after carrying out the Wick contractions for 
the graphs in fig. 2 is shown in appendix B. As a main qualitative difference with respect 
to the continuum computation, we note that the tilted Wilson lines in fig. 3 do not cancel 
against each other even at distance r± = 0. Rather, we obtain an "intercept" which we 
denote by 

2( g^\ 3 



i{v)^'-^r'-^^^, k^vp, p^yzkf, (3.9) 



where J^. and ki are defined in eq. (B.l). Then we expect eq. (3.5) to be replaced through 

1 



v}^\oo,l,r^)c^-i9'TcJlil)+ f (1 - cos V±) (7^-7-2 ^,2^ 2 

L Jkx \kfi + ki kfi + ki+ mi 

where k^ = {ky,kz), and the Debye mass parameter reads [23, 24, 39, 40] 



(3.10) 



ml = 2g^TN,— , i: = T^[^]T^[—]^^^-^ . (3.11) 

^ ^ 47ro ' ^24^ ^24^ 48^2 V ; 

As discussed in appendix B, the r_|_-dependent part of eq. (3.10) is an approximation, but is 
expected to be valid for r_|_ S> a. In any case eq. (3.10) illustrates the general feature that, 
apart from the scale of the lattice spacing, the potential can have non-trivial structure only 
at two distance scales, namely l/g'^T and l/m^. 
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Figure 4: Left: The intercept from eq. (3.9), in units of eq. (3.15), compared with lattice data. Right: 
The potential from eq. (3.10) ("LO" ); with an integrand taken from ref. [13] ("NLO" ); and after adding 
the long-distance asymptotics [15] ("NLO-I- asymptotics"). Vertical lines indicate distances beyond 
which perturbation theory is unreliable. The asymptotic behaviour sets in at f > 2rQg'^T/{2Nc) ~ 0.73. 
Axis ranges have been chosen to agree with fig. 7(left) in which lattice data is shown. 



3.3. Beyond perturbation theory 

Let us extract lessons from above for what we may expect to see in the simulations: 

• At "short" distances, r±_ <^ l/m^, ImV^i should start off with a non-zero intercept, 
given by eq. (3.9) for large Pq. 

• At "intermediate" distances, r± ~ ^/m^, the potential ImV^i shows a non-trivial struc- 
ture which is relevant for jet quenching. This structure cannot be studied quantitatively 
with the approach of the present paper, given that within CLGT the Debye scale is com- 
pletely determined by lattice artifacts, cf. eq. (3.11). On the qualitative level, however, 
we expect large corrections to the leading-order expression in eq. (3.10) [13]. 

• At "long" distances, r_|_ ^ l/m^, the phenomena related to the Debye scale are expo- 
nentially screened, and the physics is dominated by the colour-magnetic scale g'^T/ir. 
More precisely, in continuum the imaginary part of the light-cone potential corresponds 
to the static potential of three-dimensional pure Yang- Mills theory [15], which for 
r± > Tr/ig'^T) evaluates to |ImFci| ^ 0.553{g'^Tfrj_ for A^c = 3 [41]. When summed 
together with the NLO result from ref. [13], which already includes a part of the linear 
term, the appropriate correction reads 6\ Im V"ci| ^ [0.553 — 7/{l67r)]{g'^T)'^r± [15]. 



In order to observe the features mentioned in the data, it is helpful to change units. Suppose 
that we use /?g from eq. (3.7) in order to convert lattice units to physical units. Then we can 
express distances and the potential as 



r^g^T r_L 



2N, a/3G ' 

ImFci _ 1 2iVcIniFci _ t^G^maVd 
16 ^ 16 'g^ ~ 16 ' 

where the factor 16 is a convention. In these units, the Debye scale corresponds to 

1 _ g^T I vr 



friE 2iVc"iE V I^G^^c 

the zero-distance intercept from eq. (3.10) amounts to 



(3.12) 
(3.13) 

(3.14) 



hm \J^I^^^J^,C,m^m, (3.15) 

f<l/mE 16 16 2 

whereas the long-distance asymptotics reads 

lim * ™ "^^ Rs — 0.553(2iVc)^f ^ 1.2f . (3.16) 

f>l/mE 16 16 

The various features together with the effect of NLO corrections [13] are illustrated in fig. 4. 
The scale I/ttie defines the point beyond which perturbation theory is no longer to be trusted, 
and the result may eventually (for f>l) go over into the "asymptotics" curve reflecting 
non-perturbative colour-magnetic dynamics. The non-perturbative contribution of the scales 
k± ~ g'^T/iT to the so-called jet quenching parameter, g, is however determined by distances 
just above r±_ ~ l/m^^ rather than by the long-distance asymptotics [15]. 



4. Numerical implementation 

The simulations in the CLGT formalism are carried out along the lines of ref. [25], in which 
the case u = was considered. With the choice of temporal gauge t/o(x, t) = 1 the classical 
equations of motion for the variables of eq. (3.8) read 



adtUi{yL,t) = 


= i{2N,)-2Si{-K,t)Ui{-K,t) , 


(4.1 


a5t^i'(x,t) = 


/ 2 \5 

= - — ImTr 


T'u,{^,t)y^si{^,t) 


(4.2 



where Sij denotes a staple. These diff'erential equations are solved on a hypercubic three- 
dimensional spatial lattice using the Euler forward finite-difference scheme with temporal 
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Figure 5: Left: Time dependence of the classical Wilson loop for w = 1 at different distances r. The 
common fitting range for the determination of Im Vd is denoted by the shaded region. Right: Effective 
mass plots for three selected f. The fitting range is determined such that a satisfactory signal-to-noise 
ratio is obtained, however systematic errors could be substantial and the results obtained should be 
thought of as upper bounds as usual. 

lattice spacing at = -^. As initial conditions, we deploy field configurations thermalized 
according to refs. [17, 18, 19, 25] with the Hamiltonian of eq. (3.8) and the appropriate 
projection to the hypersurface respecting the Gauss law. 

To obtain the potential of eq. (3.4), we measure the discretized Wilson loop (defined like in 
eq. (2.3) but in Minkowski signature) in real time for several different transverse separations f. 
The tilting away from the temporal axis with velocity v is implemented as indicated in fig. 3. 
Fitting the time evolution of these purely real quantities (left panel of fig. 5) with an expo- 
nential allows us to read off Im V^i from the exponent. For a rough estimate of limi_!.oo Im V^i, 
we identify a common fitting range for all values of f, in which the asymptotic exponential 
falloff appears to have stabilized, while at the same time the statistical noise due to a finite 
number of measurements is still relatively small (right panel of fig. 5). These requirements 
are hard to satisfy for large f and /3g and, as can also be deduced from fig. 5(right) by bare 
eye, the procedure adopted is likely to lead to an overestimate of limt_>oolni V^ci- 

The effect of a finite volume on the determination of Im V^i is shown in fig. 6 for v = 1 and 
V = 2. Higher velocities lead to a faster exponential damping of the Wilson loop, hence the 
region for an exponential fit shrinks and leads to a more noisy signal as shown in the right 
panel. We find that to go to f ~ 0.4 a lattice extent of at least N > 1.5/3q is necessary. 

Once an adequate lattice extent and a usable fitting range have been established, we proceed 
to measure the velocity dependence of ImV^i and its intercept. For f^l/rh^ perturbation 
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Figure 6: Volume dependence of the imaginary part of the potential as obtained from lattice simu- 
lations with /3q = 64, z; = 1 (left) and v — 2 (right), and a fitting procedure as described in fig. 5. 

theory becomes more accurate at larger /3q, and indeed the lattice results approach the 
perturbative ones for the intercept at all velocities, cf. fig. 4(left). 

For f > l/m-E, in contrast, perturbation theory need not be accurate. It is perhaps surpris- 
ing then how well the "NLO+asymptotics" result works for moderate /3g, cf. figs. 4(right) and 
7(left), even though discrepancies remain at the smallest and largest /3g's- For the smallest 
/3g's this may be due to the fact that the Debye scale is larger [cf. fig. 4(right)] and therefore 
the asymptotics is approached at larger separations. For the largest /3q we reiterate that it is 
difficult to reach the regime t/a > /Sq needed for extracting the correct asymptotics (cf. fig. 5), 
so that the remaining discrepancy is probably due to systematic errors. Comparisons with 
perturbation theory need to be refined with other methods in the Euclidean domain [13], in 
which the Debye scale is free from lattice artifacts, so that the infinite volume and continuum 
limits can be systematically taken. 



5. Conclusions and outlook 

The purpose of this exploratory study has been to probe the contribution that highly occupied 
classical gauge fields make to the thermal expectation value of a light-cone Wilson loop. We 
have observed that discrepancies to the leading-order expression set in already at relatively 
short transverse distances, and lead to a larger magnitude of the imaginary part of the 
potential (stronger interactions) than predicted by leading-order perturbation theory (cf. 
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Figure 7: Left: The /3g dependence of the imaginary part of the potential for w = 1 at a fixed physical 
volume N/Pg = 1.5. Right: Velocity dependence for [3a = 64 and N//3g = 1.5. Only statistical errors 
are shown; systematic errors are expected to be larger than the statistical ones particularly for large 
/3g and imply that results should in general be interpreted as upper bounds for the true answers. 

fig. 4(right) vs. fig. 7(left)). This is in qualitative agreement with the NLO computation 
of ref. [13] and with the long-distance asymptotics as analyzed in ref. [15]. Quantitative 
comparisons are hard because of discretization artifacts inherent to the CLGT framework. 

In addition, we have noted that crossing the light cone does not change the structure of 
the potential in any qualitative way (cf. fig. 7(right)). This poses well for the proposal of 
ref. [13] according to which the potential could be measured within a purely static dimen- 
sionally reduced effective field theory [42, 43]. Unlike classical lattice gauge theory, that 
framework is (super)renormalizable, so that divergences and discretization artifacts can be 
handled through local counterterms and analytic computations, and the genuine continuum 
physics of the momentum scale k± ~ mg can be disentangled. (It is useful to stress again that 
asymptotically large values of r± need not be studied [15].) Thereby the existence of large 
infrared effects contributing to jet quenching can possibly be confirmed, perhaps leading to 
a QCD-based explanation for the experimentally observed efficient jet quenching in current 
heavy ion collision experiments at the LHC. 

We would finally like to pose the question of whether the observable of eq. (2.3) can also 
be addressed with direct four-dimensional lattice simulations. One lesson from our study is 
that discretizing the tilted Wilson lines (cf. fig. 3) is inconvenient. It might rather be sensible 
to boost the ensemble by making use of shifted boundary conditions [34], and measure the 
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Wilson loop always along the time-like lattice direction. 

Of course, measuring eq. (2.3) is not enough, but subsequently analytic continuations are 
needed for extracting the proper real-time physics. In fact there are two separate analytic 
continuations here: t ^ it as well as v^ -^ —iv. The former is conventionally implemented 
by going through frequency space, i.e. estimating the spectral function corresponding to the 
Euclidean correlator; from the spectral function, any time ordering can be recovered. For 
fE = 0, a determination of the spectral function has been attempted [27, 44, 45], and even 
though systematic uncertainties remain difficult to quantify, the challenge should not be much 
harder in the presence of v^ ^ 0. Note that the quantity of interest here corresponds to the 
imaginary part of the real-time potential, cf. eqs. (1.1), (1-2). 

As far as the analytic continuation of the velocity is concerned, one of the methods used in 
studies of QCD with a baryonic chemical potential might turn out to be helpful. For instance, 
one could first carry out simulations with v^', fit the results to a Taylor series; and subsequently 
carry out an analytic continuation. Although in a mathematical sense a singularity cannot be 
excluded as w — )■ 1, we have not observed any drastic changes in the infrared dynamics of the 
system in this limit (the discretization-specific wobbles around v ~ 0.3 in fig. 4(left) are not 
expected to be present if a boosted ensemble is simulated). Therefore it is conceivable that 
such a procedure could yield at least qualitative results against which dimensionally reduced 
simulations, carried out on the space-like side of the light cone, can be compared. 
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Appendix A. Leading-order perturbative computation in continuum 

We compute the graphs in fig. 2 with the propagator of eq. (3.1), first in Euclidean space-time. 
Carrying out Wick contractions, it can be checked that any gauge parameter dependence 
cancels. Inserting eq. (3.3) for 1/{K'^ -1-11^.^,), the remaining expression reads 

cf(r,t;E,rx) = 1, (A.l) 



xiPE(fco>k) 



1 f^^k„ 



(A;„ + k-VE)2V fc2 
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+ PTikQ,'k) 



+ 



k-rx)2 {kn + k-v^y^ 



fc2 



(A.2) 



Here /c„ = 27rn//3, with n G Z, are the Matsubara frequencies. The apparent poles of eq. (A.2) 
at /c„ + k • Ve = are regulated by the zeros of the numerator. 

The Matsubara sums can be carried out by partial fractioning the dependence on /c„, and 
then making use of 



E 



^XHjiT 



kn ikr, 



^0 






tfifiT 



Kf\ Ifxifi 



n, 



,(A:o)e(^~^)'=o , < T < /3 , (A.3) 



where n^{kQ) = l/(e^ « ~ !)• Iii order to simplify the expressions we also take the classical 
limit right away; recalling eq. (2.4) and setting /i — )■ 0, the results then become 






1 



(A.4) 



g'TC, {cosk-r^-l) -° 



xJ.pE{kQ,k) 
+ PT(feo,k) 



1 



(fcg + ik • Ve)2 
r2 1 



/^2 
1 _ _0 



A;2 



(k • r_L)2 (fcp + ik • Vj 



^;2 , M 



(A.5) 



It can be observed that in the classical limit, the Matsubara sum amounts effectively to 
replacing kn through —ik^. (It would certainly be possible to keep n^{ko) in an exact form, 
cf. ref. [46] for v = 0, however only the Bose-enhanced classical term 71,3(^0) ~ T/{hk()) is 
expected to contribute to the large-t limit to be taken presently.) 

Wick rotation is carried out through r = it, v^ = —iv, and the potential is extracted from 



.(2), 



^(2) 



(0), 



idtClii{it,-iv,r^) = V}i 'it,v,r±) C%^i{it,-iv,r±) . 



(A.6) 



We obtain 



V!i'\t,v,r^) 



g^TCp I (cosk-rx-1) 



°° dkr. e«(^()+k-v)t _ g-i(fco+k-v)t 



x<pE{ko,k) 
+ pT{ko,k) 



1 



1 _ ^ 
fco + k • V V A;2 



r^(/co +k- v) 



(k-rx)2 
Subsequently the large-time limit follows from 

gi(fcg+k-v)t _ g-i(fc(j+k-v)i 



+ 



fcn + k • V 



A-2 
,2 _ zo 

fc2 



lim 

t— 5-00 



/cq + k • V 



27ri 6{kQ + k • v) 



(A.7) 



{Al 
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Carrying out the integral over kn and setting also w = 1, so that k • v — )■ feii, leads to 



vP{^,l,r±) = -ig^TC^ I (1-coskx-rx) 



J-oo TT 1 k\\ k\\ ]k\ + kl 

Here we substituted k\\ — t- —/cm for simplicity. This potential is purely imaginary and, accord- 
ing to eq. (A. 6), corresponds to an exponential decay of the light-cone Wilson loop at large 
Minkowskian times, as anticipated by eq. (1.1). 

The next step is to perform the integral over k\\. This is possible by re-expressing the 
spectral function as a discontinuity of the retarded correlator across the real axis, 

p(t, k) = g°'^» + '°"''->^g'"'"-'°"-''' , (A.IO) 

and by then carrying out the contour integral. In the literature the procedure is known as a 
light-cone sum rule [14] (see also appendix A of ref. [13]), and yields 

dk\\ p{k\\,'k) k'^. , , , , 

„j vr fell k'l_ + ki 

The retarded propagator is, in turn, the analytic continuation of the Euclidean one. 

Recalling finally that the self-energy 11^ of eq. (3.1) vanishes at zero frequency, whereas 
IIj, equals the Debye mass parameter, m\, we recover eq. (3.5). 



Appendix B. Leading-order perturbative computation on a lattice 

If the computation of appendix A is repeated in lattice regularization, then the expressions 
become a lot more complicated. For instance, employing the notation 

fc.^-sm(-), A, = cos(-), j^ = j_^^J^,^ (B.l) 

and making use of Feynman rules derived from eq. (2.5), the observable of eq. (A. 2) can 
formally be expressed as (k = (k_|_, /cy), k^ = (fcy, kz)) 

Ck\T,v^,r^) = ^/ cosV±-l / -^22 7 T 



x{p^{kQ,\<i){k\\ 



1 ^^ + 1 



(^)2Vfc^ k^ 
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+ /?T(A:o,k) 



1 1 



+ />,(A;o,k)(l-M^|. (B.2) 

If we recall, however, that after the Matsubara sum and the classical limit, A;„ gets essentially 
replaced by —iko, and that for non-zero distances and large times the contribution emerges 
from k± <?72e and ko,k\\ "^g^T/ir (cf. eq. (A. 11)), then the lattice four-momenta can to a 
good approximation be replaced by their continuum limits, 

k/j, ^ kfj, , fc^ — >■ 1 . l-D-oj 

Then the first structure of eq. (B.2) goes over into eq. (A. 2). In contrast, the second structure, 
which is linear in r and independent of r^, originates from self-energy corrections of the tilted 
Wilson lines and is specific to lattice regularization. Since this short-distance contribution 
arises from "hard" scales, there is no need for resummation; we can replace the spectral 
representations by free propagators, 

d/CoPT(E)(^) ^ 1 ^g_^^ 



, TT ko- ikn kl + A;2 

Subsequently the Matsubara sum, classical limit, and analytic continuation are taken as 
usual, which ultimately leads to the intercept of eq. (3.9). 
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